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ABSTRACT 



Offshore jets or filaments are a well observed phenomenon in the California 
Current System. Their surface and subsurface structures as well as their spatial and 
temporal variability are well documented. The methods by which these filaments are 
formed is not as well understood as their structural form. The most common 
generation theories are: (1) variations in wind stress coupled with topographic 
irregularities, (2) dynamic instability and (3) geostrophic turbulence. In this study an 
attempt is made to identify the factor(s) responsible for filament formation south of 
Cape Mendocino by numerically simulating the California Current System using a two- 
layer, non-linear, primitive equation model. It is shown that baroclinic instability is 
the primary method by which filaments are formed in this region and that other factors 
such as barotropic instability, friction and non-linearities can alter the characteristics of 
the fastest growing baroclinically unstable wave. 
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I. INTRODUCTION 



A. THE CALIFORNIA CURRENT SYSTEM 

The California Current System, CCS, has been extensively studied by I lickey 
(1979). The CCS consists of four separate currents that are spatially and seasonally 
distributed off the West Coast of the United States. They are: (1) the California 
Current, a southward flowing surface jet; (2) the California Undercurrent, a northward 
flowing sub-surface jet; (3) the Davidson Current, a poleward flowing surface jet north 
of Point Conception; and (4) the Southern California Countercurrent, a northward 
flowing surface current south of Point Conception in the California Bight (Hickey, 
1979). Figure 1.1 shows the average seasonal configurations from Ikeda and fhnery 
(1984) for the CCS. Typical surface jet velocities arc on the order of 20 to 40 cm,'s. 
The California Undercurrent has a typical strength of about 5-10 cm s. For the 
profiles on the left, where the California Undercurrent is present, the surface jet depth 
is about 150 m while for the profiles without a subsurface jet, the surface jet extends 
down to about 300 m. 




Figure 1.1 Seasonal Configurations of the California Current 
System (Ikeda and Emery, 1984). 



B. FILAMENTS IN THE CALIFORNIA CURRENT SYSTEM 

The presence of offshore jets or filaments in the CCS is well documented. 
Satellite radiometry and color scanning along with intensive in situ data such as that 
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obtained during the Coastal Ocean Dynamics Experiment (CODE) have provided large 
data bases for these features. 

From an observational standpoint, filaments are well documented. Definitive 
studies by Flament et al. (1985), Kosro (1987) and Kosro and Huyer (1986), for 
example, have greatly expanded our knowledge of filament surface and subsurface 
structures along with their spatial and temporal variability. 

Filaments consist of strong baroclinic jets flowing offshore at peak velocities of 
0.5 m s or greater (Rienecker et al., 1985), with narrow return flows to the south 
(Flament et al., 1985). The offshore extent of filaments is typically on the order of 200 
to 300 km with maxima of around 500 km (Brink and Hartwig, 1985). Filaments 
advect cold upwelled coastal waters away from the continental shelf 

Filament widths are on the order of 50 km (Flament et al., 1985). The 
wavelength associated with their occurrence at numerous points along the coast is 
roughly 100 to 400 km for the region to the north of Cape Mendocino (Ikeda and 
Emery, 1984). Formation time scales are generally on the order of weeks with lifetimes 
on the order of a month. Figure 1.2 shows the surface velocity field derived from 
satellite feature tracking for a filament studied by Flament et al. (1985). 

Filaments are seasonal in nature. Few have been observed during late winter 
through early spring. Their first appearance is off the Baja Peninsula in March. 
Between April and May they occur betw^een Point Conception and Cape Mendocino 
and by late summer they reach their maximum northw'ard extent off the Oregon coast. 
Generally, filaments have been more observed south of Cape Mendocino. (Brink and 
Hartwig, 1985) 

C. FIL.4MENT FORMATION 

The precise formation process for filaments is not as well understood as their 
characteristics once formed. Numerous theories abound with the most prevalent 
causative forces being: (1) coastally linked processes associated with variations in wind 
stress coupled with irregularities in the alongshore topography, (2) dynamical 
instabilities of the CCS, and (3) geostrophic mesoscale eddy field interactions with the 
CCS. It is possible and indeed probable that more than one theory may be responsible 
for filament formation. 

Ikeda and Emery (1984) demonstrated by a combination of observational data 
and numerical experiments that coastal irregularities off the coast of Northern 
California and Oregon are strongly correlated with filament formation. They also 
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Figure 1.2 Surface Velocity Field of a Filament (Flament et ai, 1985). 

showed that baroclinic instability w'as probably the primary causative force north of 
Cape Mendocino. Batteen (1987) demonstrated that baroclinic instability could occur 
in the CCS due to imposing an unstable jet profile similar to the early summer profile 
shown in Figure l.l. Chen (1986) showed that baroclinic instability coupled with 
topographic variations are important to the destabilization of the CCS in the vicinity 
of Cape Mendocino. Mysak (1977) also demonstrated the importance of baroclinic 
instability in a numerical study of the CCS off Vancouver Island. Mooers and 
Robinson (1984) documented an eddy pair off Northern California, while Owen (1980) 
and Huyer ei al. (1984) observed strong eddy-like features in the northeastern Pacific 
close to the continental slope. 

D. PURPOSE OF THIS STUDY 

The objective of this study is to determine the causative factor(s) involved with 
filament formation along the coast of California south of Cape Mendocino. A two- 
layer, non-linear, primitive equation model with topography is used to simulate the 
dynamics of the CCS. The results strongly point to baroclinic instability as the 
primary generation mechanism for filaments in this region. 
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The specifics of the model and the various numerical experiments are covered in 
Chapter II. The linear baroclinic instability problem is formulated in Chapter III and 
the results applied to the data from the numerical experiments to verify the claim of 
baroclinic instability to the CCS. In Chapter IV, the results are discussed and 
conclusions presented. 
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II. NUMERICAL SIMULATION OF THE CALIFORNIA CURRENT 



A. THE NUMERICAL MODEL 
1. Model Equations 

Simulation of the CCS is performed using a two-layer, primitive equation, 
semi-implicit numerical scheme initially developed by Hurlburt (1974). This scheme 
has been employed in numerous ocean mesoscalc circulation studies (e.g., Hurlburt and 
Thompson, 1980, 1982; Smith and O'Brien, 1983; and Smith, 1986). Linear test cases 
have been run for comparison with linear analytic solutions to show that this approach 
is valid (Smith and Reid, 1982). Motion in each layer is governed by a momentum 
equation (2.1) and a continuity equation (2.2). 

( V*V.+Vj_*V)Vj_ + kxfV^ = -h^VP^ + ,2 1) 

dt ^ ' 



vn- 

+ V*Vj_ = 0 (2.2) 

Where i = layer index (i= I upper, i=2 lower) 

Variables and notation are defined in the Appendix. The Iluid is assumed to be 
hydrostatic and Boussinesq. The density in each layer is constant. 1 lie ellccts of 
winds, tides, thermodynamics and ihermohaline mixing are not included. 

2. Model Domain 

A rectangular region (1100 km x gOO km) was divided into 10 km by 10 km 
squares to form the grid for the numerical model finite differencing. By rotating the 
grid through a desired angle, the model can be used to simulate the Ilow within the 
CCS as indicated by Figure 2.1. 

Topography is included by approximating the major features apparent in the 
Mendocino Escarpment region. Figure 2.2 shows the digitized topography in the 
vicinity of the Mendocino Escarpment (Chen, 1986) and Figure 2.3 shows the 
topography used in the model. Principal features of the topography used in the model 
are: (1) a shelf region of 500 m depth extending out 50 km from the coastline which is 



13 




Figure 2.1 Geographical Positioning of Model Domain. 

Solid box represents a basin rotation (y) of 90“, dashed box for rotation of 120“. 
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a uniform vertical wall; (2) a shelf slope region extending from 50 km to ~ 100 km with 
an exponentially decreasing slope; (3) a modeled Mendocino Escarpment (without a 
ridge) which separates a 3000 m depth region to the north from 4500 m water depth to 
the south. The topography rotates with the model domain. 

3. Boundary Conditions 

Using the orientation of Figure 2.1, no-slip boundary’ conditions are set on the 
eastern and western boundaries. The northern boundary condition consists of a 
prescribed inflow in the upper layer. For simulated lower layer flows to the south or 
north the lower layer boundary consists of either a prescribed inflow or a radiation 
condition (Camerlengo and O'Brien, 1980) respectively. The southern boundary 
condition is set up just opposite the northern. The upper layer always has a radiating 
condition and the lower is radiating or prescribed for southward and northward lower 
layer flows respectively. 

Obviously, these boundary conditions do not reflect totally the situation in the 
CCS. By putting a closed boundary to the west, there will inevitably be shear 
produced vorticity that is not realistic. However, it is believed that the east-west 
dimension of the domain is large enough so that the vorticity production at the western 
boundaiy’ will not affect the dynamics of the CCS which is located near the eastern 
boundary. 

B. NUMERICAL EXPERIMENTS 
1. The Reference State 

In order to filter through the many possible generating mechanisms for 
filaments, variation of selected parameters is necessary. The problem is twofold: (1) 
selective variation of parameters may cause model deviations from the CCS 
configuration such that extrapolation of results from the model to the CCS is 
impossible. The results would still be important from a dynamical standpoint but 
without the connection with the 'real world' their usefulness is limited; (2) isolation of 
specific model responses by' changes in model parameter(s) can be inaccurate due to 
complex cross correlations between varied and non-varied parameters. 

In order to assess the model output variations resulting from parameter 
changes, a reference state is established. Instead of a simple configuration such as flat 
bottom barotropic flow, the reference state is selected to be as complex as possible 
involving such things as baroclinicity, topography and the presence of a pre-existing 
eddy (all suggested to be dynamically important in Chapter I). As mentioned before. 
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Figure 2.2 Actual Topography in Vicinity of Mendocino Escarpment (Chen, 1986). 
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Figure 2.3 Smoothed Topography used in Numerical Simulations. 
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applicability of the results to the CCS is vital and use of a reference state as near as 
possible to the CCS configuration keeps the study focus on CCS dynamics rather than 
general dynamics. 

The reference state consists of a two-layer geostrophically balanced, Gaussian 
jet of the form, 



hi(y) = .Ai[l-exp(-y2/2L2)] i= 1,2 (2.3) 

where L is the e-folding width scale of the jet and the y-axis oriented east-west as 
shown in Figure 2.1. The flow strength is determined by the amplitude coefficient A- 
which determines the magnitude of the surface and interface height anomalies through 
h-(y). Additionally, a two-layer geostrophically balanced, axisymmetric, Gaussian eddy 
is included. 



h.(x,y) = Ai[l-exp(-R2/2L2)] i= 1.2 (2.4) 

where R^ = (x^ + y^) and represents the radial distribution of the eddy. 

Mean upper layer thickness is set at 150 m (H^) which corresponds to a first 
internal Rossby radius, R^^, of ~ 18 km. This is comparable to the ~25 km R^ for the 
CCS (Brink and Hartwig, 1985). H 2 is dependent on the topography depth, D(.x,y), by, 

H 2 (x,y) = 4500 - Hj - D(x,y) (2.5) 

The jet e-folding scale, L, is set at 25 km. The jet is fixed parallel to the coast over the 
slope region and the eddy positioned offshore from the jet. The spatial relationships 
between the jet, eddy and topography are indicated in Figure 2.4. In Figure 2.4 and in 
all subsequent figures of model output fields, tick marks along the x and y axis 
represent a distance of 275 km and 200 km respectively. The lines of constant surface 
height anomaly can be thought of as streamlines representing the barotropic part of 
the current. The interface anomaly is sloped opposite the surface with an amplitude 
such that its baroclinic component in the opposite direction sets up the selected flow 
magnitude in the lower layer. 

The initial jet and eddy strength is set at 0.4 m/s (southward for jet and 
anticyclonic for the eddy) in the upper layer and 0.01 m/s in the same direction in the 
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Figure 2.4 Reference State Flows and Topography. 



Reference state initial conditions of (a) surface height anomaly(ctv>) and (h) height pi 
bottom above 4500 m (m). Contour intervals are'' (a) 2.5 cm and (b) 250 m. lick 
marks represent 275 km m the .x direction and 200 km in the y direction. 
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lower layer. The horizontal Laplacian frictional coefficient, Aj^, is set at 500 m^/s 
which is roughly in the middle of the range of (100-1000 m“/s) typically used in 
numerical studies. A summary of model parameters for the reference state is shown in 
Table 1. 



TABLE 1 

REFERENCE STATE MODEL PARAMETERS 
PAR.AMETER SYMBOL VALUE 


E-W Basin Extent 


Lv 


800 km 


N'-S Basin Extent 


Lx 


1 100 km 


Initial Upper Layer Thickness 


H, 


150 m 


Initial Lower Layer Thickness 


H,(x,y) 


Variable 


Maximum Basin Depth (F 


^1 ^2 Wax 


4500 m 


Depth of Bottom 


D(x,y) 


Variable 


Coriolis Parameter (@ 40* N lat.) 




0.94 X lO*"^ s‘^ 


df'dy 


Po 


2.0 X 10*11 j^-lj-l 


Gravitational Acceleration 


g 


9.8 m s'- 


Reduced Gravity 


g' 


2.0 X 10*2 ms'2 


Horizontal Eddy Viscosity Coef. 


^h 


5.0 X iq 2 s'l 


Time Step 


At 


3600 s 


Grid Increment in x-direction 


Ax 


10 km 


Grid Increment in y-direction 


Ay 


10 km 


Jet Maximum, Upper Layer 


Wj 


0.4 m s'l (southward) 


Jet Maximum, Low’er Layer 


V, 


0.01 m s'l (southward) 


Eddy Maximum, Upper Layer 


V,e 


0.4 m s'l (anticyclonic) 


Eddy Maximum, Lower Layer 


V,e 


0.4 m s'l (anticyclonic) 


Jet Position Offshore 


Lj 


50 km (from east wall) 


Eddy Position Offshore 


Le 


200 km (from east wall) 


Basin Rotation Angle 


Y 


90* (from east) 



Figures 2.5 and 2.6 show the reference state model output fields for height 
anomalies and upper layer velocities respectively. For this and all subsequent model 
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output field plots, surface height anomaly is contoured in cm, interface height anomaly 
is contoured in m and the velocities are contoured in cm/s. Based on the directional 
characteristics and the magnitudes of the flows, it is apparent that some phenomenon 
very’ similar to filament formation is occurring. Of particular note is the velocity 
structure of the offshore flow versus the inshore flow. The offshore flow has strong 
unidirectional flow while the inshore flow is slower and has more directional variability. 
The maximum offshore velocity by day 30 is about 50 cm/s and the inshore flow 
maximum about 40 cm/s. This type of structure and magnitude of flow correlates well 
with observed data in filament studies (Flament et ai, 1985; Kosro, 1987). 

There are also two important signatures of baroclinic instability in the output 
fields; (1) wave growth in the initially nearly quiescent lower layer is indicative of 
vertical energy transfer from the upper to lower layer, and (2) there is a 180*’ phase 
shift between the waves in the upper layer and the waves in the lower layer. 

2. Variation of Parameters 

Structures shown in Figure 2.5 are consistent with observations of the form 
and magnitude of filaments. In order to isolate potential causes for filaments, the 
assumptions that the wave structures in Figure 2.5 represent filament formation and 
growth are made. Since the reference state is as complex as possible, the variations of 
the numerical simulations are largely one of simplification in specific areas of the model 
to see whether filaments are altered or made non-existent. 

All of the numerical simulations are for 30 days and certain parameters from 
Table 1 are never varied. Table 2 provides a list of those parameters which are varied 
in the different experiments and what their reference state value is. 

a. Experiment No. 1 Increased to 250 m) 

Upper layer thickness was increased to 250 m to determine the model 
sensitivity to the vertical structure. Figure 2.7 shows the model output fields for 
surface and interface height anomalies. From the fields, it is apparent that the wave 
structures are significantly different than those in the reference state. They are weaker 
and seem to stop growing halfway through the sequence. This points to the 
importance of the vertical structure to the wave development process. 

b. Experiment No. 2 {Eddy Removed) 

The eddy was removed in this experiment to determine whether it was 
responsible for the wave structures by its interaction with the coastal jet. Figure 2.8 
shows the output fields for surface and interface height anomalies. The results are very 
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Figure 2.5 Model Output Fields of Height Anomalies 
for Reference State. 

Height anomaly fields for (a) surface at day 15 (cm), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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(a) 







(b) 



Figure 2.6 Model Output Fields of Upper Layer Velocity 
for Reference State. 

Velocity fields (cm/s) for (a) upper layer at day 15 and (b) upper layer at day 30. 
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TABLE 2 

VARIATION OF MODEL PARAMETERS 



PARAMETER VALUES (Notation from Table 1) 



E.\p. 4 
ref 


Hi 


D(x,y) 


''li 


''2i 


Vie 


V2e 




Le 


Y 


150 


1 


-.4 


-.01 


-.4 


-.01 


50 


200 


90 


1. 


250 


1 


-.4 


-.01 


-.4 


-.01 


50 


200 


90 


-> 


150 


1 


-.4 


-.01 


ne 


ne 


50 


200 


90 


3. 


150 


1 


-.4 


-.01 


-.4 


-.01 


150 


300 


90 


4. 


150 


0 


-.4 


-.01 


-.4 


-.01 


50 


200 


90 


5. 


150 


1 


-.4 


-.05 


-.4 


P 

1 


50 


200 


90 


6. 


150 


1 


-.4 


.01 


-.4 


.01 


50 


200 


90 


7. 


150 


1 


-.4 


-.01 


-.4 


-.01 


50 


200 


0 


8. 


150 


1 


-.4 


-.01 


-.4 


-.01 


50 


200 


180 


9. 


150 


0 


-.4 


-.01 


ne 


ne 


50 


200 


90 


10. 


150 


1 


-.4 


-.01 


ne 


ne 


50 


200 


120 


11. 


150 


1 


-.2 


-.01 


-.2 


-.01 


50 


200 


90 


12. 


150 


1 


-.1 


-.01 


-.1 


-.01 


50 


200 


90 



Notes: (1) Depths in m 

(2) Lengths in km 

(3) Velocities in m/s, + to north, • to south 

(4) Bottom depth: 1 topography, 0 flat bottom 

(5) Eddy rotation: CW for + vel., CCW for - vel.; ne= no eddy 

(6) Basin rotation angle (y) in degrees from east 



similar to the reference state. The wave structure is more uniform in this run, 
indicating that the eddy is not causative. However, it does affect the wave structure, 
c. Experiment No. 3 {Jet Moved Farther Offshore) 

The jet and the eddy were moved out farther away from the coast to see if 
wave structures were being generated as a result of boundary interactions with the 
eastern, no-slip boundary. Figure 2.9 shows the output fields for this run. In general, 
the results are very comparable to the reference state. It looks as though the full 
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(c) (d) 



Figure 2.7 Experiment # I (Hj Increased to 250 m) 
Model Output Height Anomalies. 



Heieht anomaly fields for (a) surface at dav 15 (cm), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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Figure 2.8 Experiment # 2 (Eddy Removed) 

Model Output Height Anomalies. 

Height anomaly fields for (a) surface at day 15 (cm), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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growth and extent of the waves to the west is now being hindered but their structure is 
consistent. Apparently, the boundaries can restrict but cannot generate the 
disturbances. 

d. Experiment No. 4 ( Topography Removed) 

In this case, topography was removed and a flat bottom at 4500 m depth 
utilized. There is substantial circumstantial evidence linking filaments to coastal 
topographic features such as the Mendocino Escarpment. By removing the 
topography, it is hoped that it can be determined whether bottom variations are 
causing the filament formation. Figure 2.10 shows the output fields for this run. The 
fields are nearly identical with the reference run. Apparently, topography is not 
inducing the filament formation process for this weak low’er layer flow. 

e. Experiment No. 5 {Increased Lower Layer Velocity) 

In this run the low'er layer velocity was increased from 1 cm/s to 5 cm/s 

flowing southward, which is considered to be on the order of the maximum observed 

velocity in the deep water south of Cape Mendocino (Stabeno and Smith, 1987). This 

reduces the vertical shear but increases significantly the model coupling with the 
bottom. Although the previous run tended to discount the importance of the bottom, 
it is important to see what possible effects it might have in a worst case situation. 
Figure 2.11 shows the output fields over the domain. The results are nearly identical 
w’ith the reference state with the amplitudes slightly smaller. This indicates that the 
model is relatively insensitive to southward lower layer velocity. 

/. Experiment No. 6 {Lower Layer Flow Direction Reversed) 

In this experiment, the lower layer flow w'as reversed so that it was flowing 
opposite the upper layer flow. This effectively increases the vertical velocity shear and 
examines a possible configuration of the CCS where there is a deep, weak northw’ard 
flow. Figure 2.12 shows the output fields for this run. The wave structure is similar 
to, yet noticeably stronger than, the reference. This would indicate the wave growth to 
be very sensitive to the vertical velocity shear. 

g. Experiment No. 7 {Zonal Flow Westward) 

In order to perform a sort of check on the model dynamics, the basin was 
not rotated into its CCS position. Since most instability studies are concerned with 
zonal rather than meridional flows, the analytic analysis of this type of problem is 
common and the expected dynamics are supported by observational data. With y = 0, 
the jet is w'estward and the eddy rotates as before, anticyclonically. Figure 2.13 shows 
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Figure 2.9 Experinrent # 3 (Jet Moved Farther Offshore) 

Model Output Height Anomalies. 

Height anomaly fields for (a) surface at day 15 (cm), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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Figure 2.10 Experiment # 4 (Topography Removed) 

Model Output Height Anomalies. 

Height anomaly fields for (a) surface at dav 15 (cm), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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Figure 2.11 Experiment # 5 (Increased Lower Layer Velocity) 

Model Output Height Anomalies. 

Height anomaly fields for (a) surface at dav 15 (cm), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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(c) (d) 



Figure 2.12 Experiment # 6 (Lower Layer Flow Direction Reversed) 

Model Output Height Anomalies. 

Height anomaly fields for (a) surface at dav 15 (cm), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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the output from this simulation. It is readily apparent that there are some type of 
instabilities in the flow but they are substantially reduced from the cases where there 
was meridional How. This makes sense from a linear instability standpoint and will be 
clarified in the next chapter. 

h. Experiment No. 8 {Zonal Flow Eastward) 

As with the previous experiment, the basin was not rotated into its 
standard CCS orientation. To examine the model response for a jet flowing zonally 
eastward, the basin was rotated 180°. As mentioned before, the dynamics of this type 
of configuration are well known and should help verify the model responses. Figure 
2.14 show's the output fields for this run. Of note is the fact that the flow is almost 
completely stable throughout the model run. Again, this makes sense from a 
dynamical standpoint and will also be addressed in the following chapter. 

i. Experiment No. 9 {Eddy and Topography Removed) 

In experiment # 2 the eddy was removed and in experiment # 4 
topography was eliminated. This does not totally resolve their involvement with the 
wave structure development. A requirement for baroclinic instability is a perturbation 
of the mean flow'. It is conceivable that the eddy and the topography could be causing 
the jet to be deflected so as to undergo perturbation growth by instability theory. To 
see whether this is true, both the eddy and topography were removed for this run. As 
can be seen in Figure 2.15, the removal of both the eddy and topography did not 
hinder at all the development of the wave structures. There are some slight differences 
from the reference state run, but they are so small that it would be difficult to 
speculate as to their exact connection w'ith the altered conditions of this run vice the 
reference one. 

j. Experiment No. 10 {Northwest to Southeast Flow) 

As can be seen in Figure 2.1, the coastline south of Cape Mendocino 
changes orientation from north-south to northw'est-southeast. To examine the effects 
of this orientation on the flow the basin was rotated 120° as shown in the figure by the 
dotted box. This new configuration is important to the degree of applicability the 
model results will have to the CCS. It is obvious that the coastal jet cannot continue 
straight south after Cape Mendocino and if wave structures are not observed W'ith the 
modeled jet going southeast it is doubtful whether the meridional flow results are true 
representations of the dynamics of the CCS. Figure 2.16 shows the model response to 
the basin rotation of 120°. The eddy was removed from this run so the wave structure 
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Figure 2.13 Experiment # 7 (Zonal Flow Westward) 

Model Output Height Anomalies. 

Height anomaly fields for (a) surface at day 15 (cm), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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Figure 2.14 Experiment # 8 (Zonal Flow Eastward) 

Model Output Height Anomalies. 

Height anomaly fields for (a) surface at day 15 (cm), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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(d) 



Figure 2.15 Experiment # 9 (Eddy and Topography Removed) 

Model Output Height Anomalies. 

Height anomaly fields for (a) surface at dav 15 (cm), (b) surface at day 30 (cm), (c) 
interlace at day 15 (m), and (d) interface at day 30 (m). 
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changes could be better seen. There are substantial differences from the reference but 
they are clearly in terms of the wave structure direction and wavelength rather than 
representing a totally different type of phenomenon. In the next chapter, the reasons 
for these structural differences will be examined. 

k. Experiment No. 11 ( Fj j Reduced to 0.2 mjs) 

Up until now all of the model runs have used an upper layer velocity of 0.4 
m, s. This may be representative of the CCS at certain times but by no means is the 
CCS consistent. As with the previous experiment, applicability to the real world 
demands that the model parameters be as realistic as possible in both value and range. 
The velocity is reduced to 0.2 m/s thus effectively halving the vertical shear. Figure 
2.17 shows the results of this alteration. The wave structures are definitely present in a 
form similar to the reference although their strength is considerably less. This makes 
sense from an instability standpoint as the shear and thus available potential energy of 
the system has been reduced. This wll be discussed in the next chapter. 

/. Experiment No. 12 (Fjj Reduced to 0.1 mjs) 

As a follow on to the previous experiment, the upper layer velocity is 
reduced again by fifty percent. The results are indicated in Figure 2.18. There is some 
very slight wave development but it is considerably less than either the previous case or 
the reference state. This would intuitively mean that the reduction of the upper layer 
velocity is causing the system to approach a point of stability. The fact that the 
system is becoming stable wth a shear still present will be covered also in the next 
chapter. 
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Figure 2.16 Experiment #10 (Northwest to Southeast Flow) 

Model Output Height Anomalies. 

Height anomaly fields for (a) surface at dav 15 (cm), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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Figure 2.17 Experiment #11 (Vjj Reduced to 0.2 m's) 

Model Output Height Anomalies. 

Height anomaly fields for (a) surface at dav 15 (cm), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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Figure 2.18 Experiment #12 (Vjj Reduced to 0.1 m/s) 

Model Output Height Anomalies. 

Height anomaly fields for (a) surface at dav 15 (cm), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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III. BAROCLIiNIC ANALYSIS 



A. SETTING UP THE LINEAR INSTABILITY PROBLEM 
1. Instability of a Non-Zonal Flow 

Most examinations of barodinic instability start with the assumption of a 
zonal flow. This assumption cannot be made when dealing with the CCS where the 
basic state is definitely non-zonal and in some regions completely meridional. To 
develop an analytic model for the CCS, the quasigeostrophic equations of motion will 
be employed. This represents a divergence from the model dynamics which arc 
operating under the primitive equations. To evaluate whether the quasigeostrophic 
equations are reasonable approximations, the Rossby number, Ro, can be calculated. 
Using U = 0.4 m's, f=0.94xi0‘'^ and L=R^i=l8.4 km yields Ro = 0.23 . This 
indicates that nonlinearities arc appreciable but the quasigeostrophic equations should 
be within the required accuracy. 

A non-zonal flow automatically implies the existence of an external forcing 
field (Pedlosky, 1979). In this case, it is assumed that the external forcing involved 
with the CCS is wind stress curl. For a two layer system without bottom (fiction the 
governing quasigeostrophic potential vorticity equations in non-dimensional form from 
Pedlosky (1979) arc: 
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(3.2) 



Notation is dedned in the Appendix. 

Following the derivation by Pedlosky (1979) for the phase speed of a 
barodinic wave in a non-zonal flow, consider the following basic state: 



= '*'i = VjX - Ujy 



(3.3) 
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= 'i'l = 0 -- 



(3.4) 



The parameters and are constant velocities in the upper layer in the x and y 
direction respectively. Note that the basic state is characterized by a uniform flow 
which does not vary in the cross stream direction. This will inherently eliminate 
barotropic instability from the analytic development but due to the length scales of the 
instabilities observed in the model runs (L> > Rj). barotropic instability is not likely 
to be a dominant force (Kamenkovich et ai, 19S6). There are no boundary constraints 
on the basic flow eliminating shear induced vorticity which is present in the numerical 
model. Again, it is felt that this will primarily affect barotropic instability modes. The 
wall ellccts on baroclinic instability are considered negligible. As Mysak ci ai (1981) 
pointed out. 

. . . inclusion of vertical side walls at the edges of a coastal current flowing along 
a continental slope does not significantly affect the growth rate, frequency or 
wavelength of the maximum unstable wave .... 

The real motivation for making these assumptions is that the analytic modeling of a 
non-uniform, non-zonal flow with boundary effects would be far beyond tlie scope of 
this presentation and would probably not yield significantly different results. .Model 
runs with small flows in the lower laj’er in both a north and south direction were so 
similar that the assumption of no flow' in the lower layer should be consistent with the 
numerical model. 

To examine the stability of the basic state, consider perturbations of the 
following form; 



Vi = 'i'l + (p,(x,y.t) 



(3.5) 



xii, = + (j)2(x.y,t) 



(3.6) 



Substitution into 3.1 and 3.2 and linearization yields 




(3.7) 
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0 



(3.8) 



I— H ^ (h^2^i) 



(.(|), f(p, 

- F,v,— * = 

c'K * ' c y 



Note that the wind stress curl forcing does not appear in the perturbation equations. 
As Pedlosky (1979) describes it, 

. . . the stability of basic flows which are not solutions of tlie unforced equations 
of motion may be consistently considered in the context of the unforced 
perturbation equations, without the need to consider explicitly the forces required 
to produce a basic state. 

This provides the basis for not specifying the form of the wind stress curl in 3.1. 

Using the normal mode approach, plane wave solutions of the following form 
are assumed: 



'*'n ^^n °''^P i(kx + fy-at) n = 1, 2 (3.9) 

Substitution into 3.7 and 3.8 leads to 



( cr-Ujk-v,f)( K’+F, ) + pk 

+ Fj(u,k+VjC)] - A2 Fj( a-Ujk-v,f ) = 0 



A2(a(K^+F2) + pk - F^Cu^k+VjC)] - A,F2 <t = 0 (3 n) 

where K~ = 

A coordinate system is now defined as shown in Figure 3.1. Ihc angle the 
perturbation wave vector makes with the x-axis is 0 while a represents the angle the 
shear flow makes with the x-axis. 
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Figure 3.1 Coordinate System Orientation. 

Based on the flow directions of the CCS within the numerical model domain noted in 
section 1 1. A, a can take on only one of two values, 270* or 300”. The coordinate 
system allows the following transformations and the introduction of the shear velocity, 



Pk/K = p cos 0 



(ujk + VjC)/K = VjCos(a*6) 

Where + 

The perturbation wave phase speed c is allowed to take on complex values, 

c = <t/K = Cj. ± icj 

After dividing 3.10 and 3.11 by K, use of 3.12, 3.13 and 3.14 results in. 



(3.12) 



(3.13) 



(3.14) 
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A,[(c-Vgcos(a-0))(K2^.F^) + pcosO + FjV 3 Cos( a-0 ) 1 (3.15) 

- A 2 Fj( c-VgCos( u-0 ) ) = 0 

A,(c(K2 + F 2 ) + pcosG - F 2 V 3 COs(a- 0 )] - AjF^c = 0 (3.16) 

Nontrivial solutions for Aj and A 2 exist only if the determinant of their coefficients in 
3.15 and 3.16 is equal to zero. Since the basic intent of this analytic development is 
specifically to support analysis of the model output and the CCS it is more practical at 
this point to switch to dimensional form. Using the notation of Kamenkovich et al. 
(19S6), the zero determinant condition leads to the phase speed of the perturbation 
wave (Equation 3.17). 

V3K2(K-+2A,)cos(a-6) - P( 2K-+A.+A, )cosO 

c = ! 1 

2K^(K'+Aj+A2) 

[ ( Aj+A2)-p‘cos20 + 2pV3K-‘( Aj-A2)cos0cos(a-0) (3.17) 

- K‘^V^-(4A,A,-K-‘)cos-(a-0)] 1/2 

± ■ ■■ ■ ■ * ' > 

2K“'( K'^+Aj+A-^) 



Where Aj = Fj/L^, A 2 =F 2 /L^, P = Pg and all velocities from before multiplied by 
scaling velocity U. 

2. Critical Velocity and Growth Rate 

When the radicand in 3.17 is equal to zero, the wave is at the transition point 
between stability and instability. At this point, can be solved for. This results in an 
expression for the critical shear velocity for instability, 

P K^(A,-A2)cos0 ± 2P(A,A2)1/2[ (A,-t-A2)--K'*] l/2cos 0 
" K^( 4AjA2-K'^)cos(a-0) 
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Utilizing the two basic state flow directions to match the CCS (a = 270°, 
300°) and using model values of Aj = 2. 94X10*^ m’^, A 2 = 1.02 x 10'^*^ m'^ and 
P = 2.0 X 10'^^ m'^’^ leaves a function of only K and 0. Figures 3.2 and 3.3 show 
the marginal stability curves for a of 270° and 300° respectively for various 
perturbation wave directions, 0. Note that curves are not shown for cases where the 
wave direction is perpendicular to the flow direction since by 3.18 this results in infinite 
shear velocities and no instability is possible. The form of the marginal stability is 
determined by the quadrant that 0 is in. The unstable region is shaded and represents 
the area where the shear velocity is greater than required for instability. The excess 
shear then becomes the energy source to provide a positive growth rate in those 
regions. The full implications of a meridional flow and a meridional wave direction arc 
indicated by Figures 3.2b and 3.2e.. The required shear flow goes to zero and thus any 
non-zero shear will result in instability as long as the wavenumber is to the left of the 
vertical asymptote corresponding to the condition, 

K < (4AjA2)‘/'* (3.19) 



which comes directly from Equation 3.18. 

With the shear conditions for instability defined, the next step is to find the 
magnitude of the growth rate in the unstable regions. Growth rate values will provide 
the necessary information to determine the characteristics of the fastest growing wave. 
For an unstable wave to grow, the radicand of 3.17 must be negative and the growth 
rate, Kc| becomes, 

[ K^Vg2( 4A,A2“K^)cos2(a-0) ^ 

_ - 2pV3K^(A^-A2)cos0cos(a-9) - ( A^-^A,) Vcos-0] 

2K(K‘+Aj+A2) (3.20) 



As before, model values for Aj, A 2 and P along with flow angle a reduces the growth 
rate dependence to just V^, 0 and K. Figures 3.4 and 3.5 show the growth rates for a 
= 270° and 300° respectively for velocity shears of 0.4, 0.3, 0.2 and 0.1 m/s. As 
expected, the growth rate magnitude is a strong function of the shear velocity since 
that provides the energy input into the growing wave. The maximum growth rate. 
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Figure 3.2 Marginal Stability Curves for a = 270*. 

Curves for 0 of (a) 45*, (b) 90*. (c) 135*, (d) 225*, (e) 270* and (Q 315*. 
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Figure 3.3 Marginal Stability Curves for a= 300”. 

Curves for 0 of (a) 75°, (b) 120° (c) 165°, (d) 255°, (e) 300° and (f) 345°. 
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(a) 



(b; 





(c) 



(d) 



Figure 3.4 Growth Rates for a = 270®. 

Growth rates ( x 10’^ s‘^) for velocity shears of (a) 0.4, (b) 0.3, (c) 0.2 and (d) O.I m/s. 



48 





(a) (b) 





(c) (d) 



Figure 3.5 Growth Rates for 0 = 300°. 

Growth rates ( x 10'^ s'^) for velocity shears of (a) 0.4, (b) 0.3, (c) 0.2 and (d) 0. 1 m/s. 
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which would correspond to the fastest growing, wave, occurs at different K-6 
intersections in each different shear flow case. This means the shear flow not only 
affects the magnitude of the growth rate but also the direction and wavenumber of the 
perturbation wave. As a result, any variation in the external forcing which is 
producing the shear flow would cause major alterations to the fastest growing wave. 
Another interesting point is that the growth rates for 0 = 270® are higher in all shear 
conditions than for a =300® although the 0-K plane positive growth rate areal extent 
is smaller. From the marginal stability curves, it is apparent that a meridional flow 
could have perturbations for which any non-zero shear would support instability. As a 
result, the additional claim that meridional flow yields the largest growth rate for a 
given shear can be made. 

B. APPLICATION OF ANALYTIC RESULTS TO MODEL OUTPUT 

The instability problem has now been formulated to approximate the numerically 
modeled situation. The next step is to examine the instabilities present in the model 
output fields as noted in section II.B.2 in terms of the previous section's results. If 
there are no inconsistencies, then the instabilities are most likely due to baroclinic 
instability of the basic flow. 

From section II.B.2, it is obvious that there is a similarity of the output fields 
with respect to apparent wave-like structures. Not only do the same features appear in 
most of the runs, but their shape, direction and magnitudes are also similar. Table 3 
shows the average values of relevant parameters from the model runs (where observed) 
for days 10, 20 and 30 for flow direction of a=270®. Table 4 shows the same 
parameters for flow direction of a= 300®. The wave data for both tables was taken 
from model runs where the shear flow was close to 0.4 m/s ( ± 0.01 m/s). 



TABLE 3 

WAVE PARAMETERS FROM MODEL OUTPUT FIELDS (a = 270®) 



DAY 


e (deg.) 


K(m-l) 


L (km) 


10 


215 


2.0 X 10-^ 


310.0 


20 


249 


2.6 X lo*^ 


240.0 


30 


255 


X 

o 


200.0 
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TABLE 4 

WAVE PARAMETERS FROM MODEL OUTPUT FIELDS (a=300°) 



DAY 


0 (deg.) 


K(m-l) 


L (km) 


10 


238 


2.3 X 10'^ 


275.0 


20 


260 


O 

X 

o 

i 


209.0 


30 


269 


3.2 X 10-5 


198.0 



These days were selected because they represent the appearance (day 10), the 
growth (day 20) and the maturity (day 30) of the waves. By mapping the wave data 
from Tables 3 and 4 onto instability plots from section III.A.2, an evaluation of their 
relationship to the analytically derived baroclinic wave characteristics can be made. 

•Marginal stability curves are shown in Figure 3.6 for the flow angles and wave 
angles noted in Tables 3 and 4. The wavenumbers for those angle values are indicated 
and marked by the day the data is associated with at its intersection with the Vj = 0.4 
m/s line. Of note is the fact that all of the plots show the waves within the unstable 
region for baroclinic instability. Additionally, as the waves shift in wavenumber 
towards shorter wavelengths, their shear excess increases to its maximum where the 
wave meets the asymptote corresponding to Equation 3.19. 

Growth rate plots are shown in Figure 3.7 for velocity shears of 0.4 m/s. The 
parameter positions corresponding to days 10, 20 and 30 from Tables 3 and 4 are 
indicated to show the path the waves took as they developed. From Figure 3.7 it is 
apparent that the model output waves fall in the regions of positive gromh rate but do 
not fall on or near the peak growth rate which would correspond theoretically to the 
position of the fastest growing wave. Also, the growth rate is low at day 10, reaches a 
ma.ximum at day 20 and is decreasing by day 30. It must be kept in mind that these 
growth rates represent the analytic values given the wavelength and wave direction of 
the model output waves. They are not necessarily the same as the growth rates of the 
actual waves in the model output. The basic effect of the growth rate is to increase the 
amplitude of the growing wave in the following manner: 

A = Age^^^) (3.21) 
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Figure 3.6 Marginal Stability Curves for Model Output Wave Characteristics. 

Curves for a = 270® and 0= (a) 215®, (b) 249®, (c) 255® and for a = 300® and 0= (d) 
238^(e) 260®, (0 269^ 
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Figure 3.7 Growth Rate Positions of Model Output Waves for Shear Velocity of 0.4 m/s. 

Parameter data from (a) Table 3 and (b) Table 4 ploted on growth rate ( x 10'^ s*^) 
fields from Figures 3.4. a and 3. 5. a. 

where A is the time varying amplitude and Aq is the initial amplitude. Rearranging 
Equation 3.21, 



Cn(A - Aq) = kcjt (3.22) 

yields an expression of the form, y= mx where m represents the slope (kcj). By plotting 
the natural log of the amplitude difference from the initial versus time, the growth rate 
can then be determined by the slope of the curve. Figure 3.8 shows the amplitude 
changes of the reference state waves during the model run. 

Immediately apparent is the fact that the growth rate varies with time. It is low 
at first and grows to a maximum and then decreases. The information from Figure 3.8 
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Figure 3.8 Reference State Wave Amplitude Growth. 

can be easily converted into a plot of growth rate versus time. This is shown in Figure 
3.9. The temporal variation of the growth rate is more clear and shows similar form to 
the growth rate progression observed in Figure 3.7.a. 

Experiments 7 and 8 involved zonal flows. These were not meant to simulate the 
CCS but to evaluate the model response against theory. For a zonal flow, the wave 
growth must offset one hundred percent of the stabilizing planetary vorticity gradient, 
whereas for a meridional flow the gradient's effects are zero due to the Pcos6 term. 
Figure 3.10 shows the marginal stability curves for an eastward and westward zonal jet. 
The curves are not for any specific 0 because the 0 dependence cancels out of equation 
3.18 for a zonal flow. Inamediately apparent is the fact that the easterly jet should be 
more stable than the westerly jet since it has a much higher minimum shear velocity. 
This agrees with the experimental results where the easterly jet remained stable through 
the entire run and the westerly jet developed slight instabilities. Based on the curves 
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Figure 3.9 Growth Rate Variations of Reference State Waves. 

for the zonal flow coupled with the fact that the shear velocity was 0.4 m/s, one might 
expect instability wave growth for both zonal flows similar to the reference run 
although this was not observed. The increase in the required for the zonal flows 
due to the P effect cuts down on the excess shear available for growth and apparently 
this excess is not quite great enough for wave growth. This points to frictional eflccts 
and will be discussed in the next section. 

C. FURTHER ASPECTS OF THE INSTABILITY PROBLEM 

The assumptions made to support the application of baroclinic instability theory 
to the numerical model results were necessary for any quantitative analysis to be made. 
Generally, there was very good agreement between linear instability theory and the 
model output which supports the contention that baroclinic instability is in fact the 
dominant dynamical process in the model and in the CCS. Nevertheless, there were 
discrepancies between the model simulations and the analytic predictions and these 
need to be addressed. 
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Figure 3.10 Marginal Stability Curves for Zonal Flow. 

Curves are for a= 270° and flows to (a) east and (b) west. 

1. Effects of Friction 

Without friction, the vertical phase shift of a baroclinic wave opposite the 
vertical velocity shear is an indication that cj>0 and the wave should grow (Pedlosky, 
1979). However, when frictional forces are introduced, the wave can take on the same 
phase shift opposite the velocity shear and still be neutral. Instead of supplying energy 
for wave growth, the extracted energy is spent offsetting dissipation. Thus, in order to 
grow in the presence of friction, the baroclinic wave must have a greater phase shift 
than would be the case without friction. Since there is a maximum phase shift 
associated with maximum energy extraction, the growth rate of a baroclinic wave with 
friction should be less than that for a wave without friction. 

In Figures 3.4 and 3.5 it was observed that changes in the shear velocity 
resulted in rna.ximum growth rate value and 0-K positional changes. Frictional forces 
would tend to reduce the available energy for growth which is equivalent to a reduction 
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in the shear velocity. For this reason, the numerical output waves could not be 
expected to fall on the maximum growth rate parametric positions predicted by the 
non-frictional analytic development. .-Another complicating factor is the way in which 
friction was parameterized in the numerical model. By using a Aj^V^V form, friction is 
a strong function of the velocity field. Since the velocity field in the domain changed 
drastically as the baroclinic waves developed, the frictional forcing must also have 
changed drastically. Thus, the frictional alteration of the maximum growth rate 0-K 
position would take on temporal variability. This might explain one reason the model 
output waves were observed to vary in wavenumber and direction through the model 
runs. 

2. Non-linear Instability Effects 

In linear instability theory the perturbation wave is unbounded as long as it 
meets the necessaiy conditions for instability which are defined in terms of the basic 
state. Obviously, this is not realistic as the weaves grow to finite amplitude. Non-linear 
theory explains how advection affects the wave grow’th. 

As Pedlosky (1979) pointed out, it is the basic state in the absence of the 
perturbation wave that must be used as criterion for instability vice any instantaneous 
flow profile. Variations in the flow profile as the wave develops, however, does have 
an effect on the wave. When the wave reaches a point where the instantaneous mean 
flow would not support instability, the w'ave growth rate is at its maximum and will 
start to decline. Although declining, the growth rate is still positive and the wave is 
still growing even though the instantaneous mean flow might support no 
instability.The change in growth rate is accompanied by a reduction of the phase shift 
of the wave in the vertical until growth is stopped as the phase shift reaches zero. The 
process does not stop here but continues into the phase orientation that causes wave 
decline. As with the growth, the decline will' hit a maximum and reverse and the whole 
process of growth and decay will oscillate back and forth (Pedlosky, 1979). Thus, the 
variations in the wave growth rate observed in Figure 3.9 and inferred from Figure 3.7 
could well be indicative of non-linear effects present in the model and presumably, the 
CCS. 

When considering the effects of friction, variations in the available energy to 
the wave were correlated to variations in the parametric characteristics of the 
maximum growth w'ave. Non-linear theory predicts similar variations in the available 
energy due to the changes in the mean flow brought about by the perturbation wave. 
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Hence, it is another possible causative factor for the. model output, since the model is 
non-linear, to exhibit variability in baroclinic wave wavenumber and direction over 
time. 

3. Mixed Instability Effects 

Most treatments of the instability problem are either focused on baroclinic or 
barotropic instability. In the CCS and the numerical model, flows have both 
horizontal and vertical shears. It was assumed in the analytic development of the 
baroclinic instability problem in section III.A.l that barotropic instability was not the 
dominant instability process based on the scales of the observed disturbances. 
However, it may be that barotropic instability does have a slight influence on the 
baroclinically unstable waves. 

Pedlosky (1979) showed that for a system with vertical and horizontal velocity 
shears, the characteristics of the most unstable wave can only be derived with 
exhaustive calculation given a detailed spatial distribution of velocity and potential 
vorticity. Changing the parameters of the flow can have a strong inpact on the type of 
instability present and the wavelength. It is believed that throughout the numerical 
simulations, baroclinic instability remained the dominant instability process but 
barotropic instability may have altered the characteristics of the fastest growing wave. 

4. Baroclinic Instability with Complex Wavenumbers 

In developing the baroclinic instability problem, a particular perturbation 
form was assumed. This form is common to nearly all studies of this type of instability 
and is shown by Equations 3.9 and 3.14. The crucial point is the selection of a real 
wavenumber and a complex phase speed. 

Hogg (1976) examined analytically and numerically the baroclinic instability 
process in a zonal flow where the perturbation form allowed both complex phase speed 
and complex wavenumber. In general, his results indicated a wavenumber cutoff 
between spatially growing (k->0) and temporally growing (cj>0) waves. Temporally 
growing waves were confined to long wavelengths while spatially growing w'aves w'ere 
confined to the shorter wavelengths. 

In a study of Gulf Stream meanders. Watts and Johns (1982) used statistical 
analysis of cross stream data and derived values for both temporal and spatial growth 
rates of the meanders. Contrary to Hogg (1976), Watts and Johns' (1982) results 
indicated that spatial and temporal growih rates can occur in similar wavelength 
regimes. Figure 3.11 shows the growth rate versus wavenumber relation for both 
spatial and temporal growth found by Watts and Johns (1982). 
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Figure 3.1 1 Spatial and Temporal Wave Growth of Gulf Stream Meanders 

(Watts and Johns, 1982). 

There is evidence both analytically and observationally that spatial growth 
may be involved in the baroclinically unstable wave. Although Ilogg (1976) 
concentrated on spatial growth (kj>0) and Watts and Johns (1982) observed only 
temporal and spatial gromh (cj>0, kj>0), there is nothing in the theory that prevents 
a combination of temporal growth and spatial decline (cj>0, kj<0). This might also 
explain the variation in the observed w'avelengths in the model output from long to 
short while the waves were growing in time. 
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IV. DISCUSSION AND CONCLUSIONS 



A. ANALYSIS OF RESULTS 

1. Comparison of Results with Observational Data 

As mentioned in the introduction, filament structure, spatial distributions and 
temporal variabilities are well known. In order to validate the results of this study, 
there should be general agreement with the existing information and if discrepancies 
are present there should be some theoretical basis for the differences. In Table 5 
values are shown for the range and magnitude of the observed versus the modeled 
characteristics. The value for the observed wavelength and direction is from Ikeda and 
Emery (1984) and the observed offshore extent is from Brink and Hartvvig (1985). The 
rest of the observed values are from Flament et al. (1985). 





TABLE 5 




FILA.MENT CHAR.ACTERISTICS, OBSERVED VS. MODELED 


PAR.AMETER 


OBSERVED 


MODELED 


Wavelength 


100-400 km 


198-310 km 


Wave Direction 


~270" 


215-269“ 


Off-shore Extent 


100-500 km 


100-800 km 


Maximum Off-shore Velocity 


~0.5 m/s 


0.5 m/s 


Maximum On-shore Velocity 


■^0.35 m/s 


0.4 m/s 


Growth rate (e-folding time) 


~ 1-2 weeks 


-10 d 



It would seem that the only parameter that deviated significantly was the wave 
direction. The reason for this will be discussed in the next section. Otherwise the 
agreement strongly supports the contention that the modeled phenomena are in fact 
representative of filaments in the CCS. 

2. The Filament Formation Process 

The filament formation process is quite complex. Although baroclinic 
instability is clearly the dominant mechanism, there are other forces at w'ork that 
cannot be neglected. 
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In the first phases of growth, the wave is -affected by the east-west horizontal 
velocity shear. Close to the jet, the wave front is advected strongly to the south. 
Farther offshore, the jet effect on the wave propagation decreases due to the mean 
velocity gradient. This causes the wave vector direction to be deflected from the 
seemingly more favorable southerly direction to a near westerly direction. As the wave 
gains energy from the vertical shear of the mean flow, advection of the mean flow in a 
zonal direction starts to break down the east-west horizontal velocity shear and the 
wave direction rotates towards its more optimum direction and its wavelength shortens 
to maximize the growth rate. Non-linearities become dominant as the wave reaches 
such a magnitude that the instantaneous mean flow is stable and the wave growth 
starts to decline. This is as far as the waves were modeled, but it would be expected 
that the waves would go into an oscillatory growih-decline mode until damped by 
friction. 

Obviously, this w'hole analysis is contingent on a constant wind stress forcing 
for periods on the order of months which is hardly realistic. It can be expected that 
constant forcing may persist for weeks at a time in which case the initial growth and 
development holds true. 

From an observational standpoint, the filaments observed are nearly always 
oriented perpendicular to the flow indicating a southerly wave direction. This does not 
mean that the described theory of development is necessarily flawed, it just means that 
by the time the waves are observable, they have probably already twisted around into 
their more optimum configuration. 

B. IMPORTANCE OF RESULTS 

There have been numerous studies of zonal instability processes and a few of 
meridional situations but there is a large and distinct scarcity of studies of flows in- 
betw'een. When one considers currents such as the Gulf Stream, the Kuroshio and, of 
course, the CCS which have flows that are not always east-west or north-south, it 
becomes evident that more attention should be placed on this discrepancy. 

As demonstrated, the behavior of a flow that is non-zonal is dramatically 
different from a zonal flow. In a zonal flow, the energy releasing horizontal 
temperature gradient and the stabilizing planetary vorticity gradient are always working 
against each other. In a meridional flow', energy release can take place without any 
component of P acting on it; hence, any shear above the dissipation level may produce 
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instability. For the in-between case, the wave must assume a balance between the 
maximization of energy release and the dampening effects of the planetary vorticity 
gradient. 

C. RECOMMENDATIONS FOR FUTURE STUDIES 

There are several aspects of this modeling work that could be modified to better 
simulate the CCS. .A better resolution model would help significantly to bring the 
barotropic instability process into view. This is particularly important to the filament 
situation since observational evidence has shown instabilities along the filament frontal 
region of the scale of that are claimed to be the result of barotropic instability 
(Flament et ai, 1985). .A model with more layers would give the vertical resolution 
necessaiy to properly simulate the CCS in all of its configurations, particularly its late 
summer profile of a strong southward surface jet and a poleward undercurrent all 
above a relatively quiescent deep water region. The addition of some temporal 
variability to the coastal jet strength would help to better understand wind stress 
variation effects on existing filaments and the growth of future filaments. Continued 
intensive in situ and satellite gathering of oceanographic data is vitally important to 
understanding the physical processes necessary to better reproduce the observed 
filament structures numerically. To thoroughly determine if baroclinic instability is 
occurring in the CCS. an elaborate experiment would have to be designed to resolve 
spatially and temporally the mesoscale motions. An experiment comparable to the 
POLYMODE experiment in the North Atlantic would be necessary. 

D. CONCLUSIONS 

Baroclinic instability is a prime generating force for filaments formed south of 
Cape Mendocino in the CCS. Barotropic instability along with non-linearities, 
frictional effects and possible spatial wave growth may help to modify the 
characteristics of the growing baroclinically unstable waves. 

Filament growth is strongly dependent on the vertical velocity shear, the vertical 
structure and the mean flow direction. Any attempt to model filaments must carefully 
parameterize these quantities if accurate and thus useable results are to be obtained. 
Indeed, for any numerical simulation of currents with unknown instabilities, their 
values are crucial to correct simulations. Improper understanding of their spatial and 
temporal characteristics can seriously affect the validity of the model results. 
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APPENDIX 

SYMBOLS AND NOTATION 






D(x,y) 

Di 

Fi 

fo 

g 

o' 

o 

hi 

Hi 

k 

k 

C 

K 

L 



Pi 

P2 

Ro 

U 



Ui,Vi 

Ui,Vi 



Vs 

Vsc 

x,y 

a 



Po 

P 



V 

I 

Ax, Ay 



Laplacian lateral friction coefficient 

Amplitude coefficient for normal mode solution for n^h 

Phase speed of perturbation wave 

Variable depth of topography 

Upper layer thickness 

Froude number for i^h layer 

Coriolis parameter for mean latitude 

Gravitational acceleration 

Reduced gravitational acceleration 

Instantaneous layer thickness 

Upper (i= 1) and lower (i= 2) layer mean thickness 

Unit vector in the z direction 

Perturbation wave number in the x direction 

Perturbation wave number in the y direction 

Total wave number for perturbation wave 

Horizontal scale length 

Pressure in the upper layer 

Pressure in the lower layer 

First internal Rossby radius of deformation 

Rossby number 

Scale velocity 

Velocities in the x and y directions 
Transport in the x and y directions 
Shear velocity 

Critical shear velocity for instability 
Cartesian coordinates directed N and W respectively 
Angle mean flow vector makes with x axis 
Variation of Coriolis parameter with latitude 
Scaled variation of Coriolis parameter with latitude 
Basin rotation angle from east 
Grid spatial resolution 



layer 



= g'H; 



= g(p2-Pi)'Pi 



= k^ + f 



= g(hj + h2 + d) 

= Pi-g'hi 

= [g'HiH2,'fo2(Hi + H2)lh'2 

= U'fU 
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At Model time increment 

6 Angle perturbation wave vector makes with x axis 

Dimensionalized Froude number for i*-^ layer 
Pj Density in i^^ layer 

Pg Mean density 

(7 Frequency of perturbation wave 

T Wind stress vector 

(p- Perturbation stream function for i^^ layer 

'Pj Basic state stream function for i^^ layer 
V}/- Variable stream function for i^^ layer 

V Gradient operator 

V" Laplacian operator 




L2 



= d;dx + d;dy 
= d^;dx“ + d~,'dy^ 
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